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Blazar Variability and Evolution in the GeV Regime 

S. Tsujimoto, J. Kushida, K. Nishijima, and K. Kodani 
Tokai university, Hadano city, Japan, 259-1292 

One of the most important problem of the blazar astrophysics is to understand the physical origin 
of the blazar sequence. In this study, we focus on the GeV gamma-ray variability of blazars and 
evolution perspective we search the relation between the redshift and the variability amplitude of 
blazars for each blazar subclass. We analyzed the Fermi-LAT data of the TeV blazars and the bright 
AGNs (flux > 4x10”^ cm“^s“^) selected from the 2LAG (the 2nd LAT AGN catalog) data base. 
As a result, we found a hint of the correlation between the redshift and the variability amplitude in 
the FSRQs. Furthermore the BL Lacs which have relatively lower peak frequency of the synchrotron 
radiation and relatively lower redshift, have a tendency to have a smaller variability amplitude. 


I. INTRODUCTION 


The blazar is the class of the active galactic nu¬ 
clei (AGN) which has the most number of extra- 
galactic source in the very-high-energy gamma-ray 
regime {E > 100 GeV)[l3. They are character¬ 
ized by double-peaked nonthermal emission with spec¬ 
tral energy distribution (SED) in radio to gamma- 
ray regime. Blazars include BL Lacertae objects (BL 
Lacs) and flat spectrum radio quasars (FSRQs). In 
addition BL Lacs include high-frequency peaked BL 
Lac objects (HBLs), intermediate-frequency peaked 
BL Lac objects (IBLs) and low-frequency peaked BL 
Lac objects (LBLs). In the leptonic model, the low 
and high bump of SED is explained by the synchrotron 
and the Synchrotron Self Gompton (SSG) and/or Ex¬ 
ternal Compton (EC) radiations. The estimation of 
EC in the second hump is a matter of great impor¬ 
tance in estimating external photons. 

One of the characteristic of the blazar spectra is 
blazar sequence. In 1998, Fossati et al. combined 
three complete blazar samples Q; The 2 Jy samples 
of FSRQs [1], the radio selected 1 Jy samples of BL 
Lacs0 and the X-ray^selected sample (Einstein Slew 
Survey) of BL Lacs @]. The thirty-third sources of 
selected sample were detected in high-energy gamma- 
ray regime {E > 100 MeV) by the EGRET instru¬ 
ment on-board the compton gamma ray observatory 
(CGRO). These sources were devided 5 bins based on 
the 5 GHz radio luminosity and averaged the SED 
of the each type of the blazars. The made SED 
suggested some relationship; first, the synchrotron 
peak frequency and the bolometric luminosity have 
the anti-correlation, second, the synchrotron peak 
frequency and the compton peak frequency have the 
positive-correlation, finally, the compton dominance 
(the ratio of the inverse Compton to synchrotron lumi¬ 
nosity) and the bolometric luminosity have positive- 
correlation. These correlations are known as “the 
blazar sequence” in considering the blazar physics. 
We aim to reveal the relation between the evolution 
process of AGNs and the blazar sequence based on the 
systematical study for many blazars . 

In this study, we calculated 100 AGNs (blazars) 


to find the difference of the variability amplitude in 
blazar types and evolution of the variability amplitude 
(activity). We applied the fractional variability am¬ 
plitude (Fvar) to calculate the variability amplitude 
considering the error. The Fvar is defined as Eq. o 
which was given by Vaughan et al. (2003). 


I 02 _ ^2 

= ( 1 ) 

Note that is the total variance of the light curve, 

_ _ 2 

is the mean square of flux error and F is the 

square of mean flux. Fvar error (uncertainty) is de¬ 
fined as Eq. m by Poutanen et al. (2008). 


AEvar — + 6rr((T^xs) ^Gar (2) 

Where err(cr^xs) i® defined as Eq.dS]) by Vaughan et 
al. (2003) 



( 3 ) 


Evar is often used in computation variability ampli¬ 
tude for each spectral bandQ- 

In this study, we studied the flux variation of AGNs 
with the Evar on the high-energy gamma-ray regime 
to get the variation charactor of each class (type). 


II. DATA SELECTION 


We selected AGNs from the second LAT AGN cata¬ 
log (2LAC)Q and TeVCat 121. Selection criteria were 
as follows, 


I : It was decided which subclass was belonged to 
(HBL or IBL or LBL or FSRQ). 


II : It had known redshift. 


Ill : flux > 4 X 10 ® cm ^ s ^ in 2LAG. 
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TABLE I: Analyzed AGN list UM 


Source 

name 

Redshift 

z 

Class 

(type) 

Source 

name 

Redshift 

z 

Class 

(type) 

Source 

name 

Redshift 

z 

Class 

(type) 

Messier 87 

0.0044 

FRI 

MS 1221.8+2452 

0.218 

HBL 

4G +55.17 

0.899298 

FSRQ 

NGC 1275 

0.017559 

FRI 

PKS 0301-243 

0.26 

HBL 

PKS 0823-223 

0.91 

IBL 

Mrk 421 

0.031 

HBL 

S2 0109+22 

0.265 

IBL 

PKS 0420-01 

0.916 

FSRQ 

Mrk 501 

0.034 

HBL 

lES 0414+009 

0.287 

HBL 

AO 0235+164 

0.94 

LBL 

lES 2344+514 

0.044 

HBL 

S5 0716+714 

0.31 

IBL 

S3 0218+357 

0.944 

Blazar 

Mrk 180 

0.045 

HBL 

OT 081 

0.322 

LBL 

OP 313 

0.997249 

FSRQ 

lES 1959+650 

0.048 

HBL 

lES 0502+675 

0.341 

HBL 

PKS 0454-234 

1.003 

FSRQ 

AP Lib 

0.049 

LBL 

PKS 1510-089 

0.361 

FSRQ 

4G +14.23 

1.038 

FSRQ 

lES 1727+502 

0.055 

HBL 

3C 66A 

0.41 

IBL 

PKS 2201+171 

1.076 

FSRQ 

BL Lacertae 

0.069 

IBL 

PKS 0735+17 

0.424 

LBL 

PKS 0426-380 

1.111 

LBL 

PKS 2005-489 

0.071 

HBL 

4C +21.35 

0.432 

FSRQ 

PKS B1908-201 

1.119 

FSRQ 

RGB J0152+017 

0.08 

HBL 

lES 0647+250 

0.45 

HBL 

OG 050 

1.254 

FSRQ 

lES 1741+196 

0.083 

HBL 

PG 1553+113 

0.5 

HBL 

PKS 1551+130 

1.30814 

FSRQ 

W Gomae 

0.102 

IBL 

GB 1310+487 

0.501 

FSRQ 

PKS 0244-470 

1.385 

FSRQ 

lES 1312-423 

0.105 

HBL 

PKS 2326-502 

0.518 

FSRQ 

PKS 2023-07 

1.388 

FSRQ 

VER J0521+211 

0.108 

IBL 

3G 279 

0.536 

FSRQ 

S4 1030+61 

1.40095 

FSRQ 

PKS 2155-304 

0.116 

HBL 

MG2 J071354+1934 

0.54 

FSRQ 

PKS 0402-362 

1.417 

FSRQ 

B3 2247+381 

0.1187 

HBL 

BZQ J0850-1213 

0.566 

FSRQ 

PKS 0250-225 

1.419 

FSRQ 

RGB J0710+591 

0.125 

HBL 

PKS 1424+240 

0.6035 

IBL 

PKS 1454-354 

1.424 

FSRQ 

H 1426+428 

0.129 

HBL 

4C 31.03 

0.603 

FSRQ 

B2 1520+31 

1.484 

FSRQ 

lES 1215+303 

0.13 

HBL 

PMN J2345-1555 

0.621 

FSRQ 

PKS 2052-47 

1.489 

FSRQ 

PKS 1717+177 

0.137 

LBL 

PKS 1244-255 

0.633 

FSRQ 

PKS 2227-08 

1.55999 

FSRQ 

lES 0806+524 

0.138 

HBL 

S4 1849+67 

0.657 

FSRQ 

TXS 1013+054 

1.7137 

FSRQ 

lES 0229+200 

0.14 

HBL 

4C +56.27 

0.664 

LBL 

PKS 0215+015 

1.721 

FSRQ 

IRXS J101015.9-311909 

0.142639 

HBL 

S5 1803+784 

0.68 

LBL 

MGl J123931+0443 

1.76095 

FSRQ 

TXS 1055+567 

0.14333 

IBL 

Ton 599 

0.724565 

FSRQ 

MG2 J101241+2439 

1.805 

FSRQ 

3C 273 

0.158 

FSRQ 

B2 0716+33 

0.779 

FSRQ 

4G +38.41 

1.81313 

FSRQ 

H 2356-309 

0.165 

HBL 

TXS 0106+612 

0.785 

FSRQ 

PKS 0805-07 

1.837 

FSRQ 

PKS 0829+046 

0.173777 

LBL 

PKS 1622-253 

0.786 

FSRQ 

PKS 1502+106 

1.83928 

FSRQ 

lES 1218+304 

0.182 

HBL 

B2 2234+28A 

0.795 

LBL 

4G 01+02 

2.099 

FSRQ 

lES 1101-232 

0.186 

HBL 

PKS 0440-00 

0.844 

FSRQ 

PKS 1329-049 

2.15 

FSRQ 

lES 0347-121 

0.188 

HBL 

OG 457 

0.859 

FSRQ 

S4 0917+44 

2.18879 

FSRQ 

RBS 0413 

0.19 

HBL 

3C 454.3 

0.859 

FSRQ 

PMN J1344-1723 

2.506 

FSRQ 

OX 169 

0.211 

FSRQ 

TXS 1920-211 

0.874 

FSRQ 




lES 1011+496 

0.212 

HBL 

PKS 0537-441 

0.892 

LBL 





Table U shows the analyzed AGN list which 102 
sources are included in. Source name, redshift, and 
class(type) are cited from 2LAC and TeVCat. 


The redshift of PKS 1424+240 was referred to 
Furniss(2013)@. According to the TeVCat (Tl|. the 
blazar class of S3 0218+357 (z=0.944) was not deter¬ 
mined but we used it as high redshift VHE gamma-ray 
emitter. 


III. DATA ANALYSIS 


We analysed the Fermi reprocessed pass 7 data 
between 2008 August 04 and 2014 June 09, using 
the unbinned likelihood analysis with the Fermi Sci¬ 
ence Tools package version v9r33p0 available from 
the Fermi Science Support Center (FSSC)[I^. The 
likelihood analysis was selected that the events with 
photon energies in the range of 0.1-300 GeV and 
a Region Of Interest (ROl) of 10 degrees cen- 
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I Mrk421 I 



PKS 0454-2^ 



FIG. 1: Examples of the gamma-ray light curve about 6 years in 0.1-300 GeV. Blue plots: data points with a bin size of 
30 days (bin size of last it is only about 4.4 days). Red allows: 95% C.L. upper limits. Green dotted lines: average flux 
of whole period. Left panel: Mrk 421(HBL). Right panel: PKS 0454-234(FSRQ). 


tered at the position of Table H] sources. We used 
“SOURCE” class (“evclass = 2”) including both front 
and back events, because the “SOURCE” class is rec¬ 
ommended for off-plane point source analysis by the 
likelihood analysis [l^. We excluded events with 
zenith angles larger than 100 degrees and time in¬ 
tervals when the rocking angle was larger than 52 
degrees. The set of the instrument response func¬ 
tions of “P7REP.SOURCE_V15” was applied. Mod¬ 
els which were used in this study include the isotropic 
diffuse background (iso_source_v05.txt [1^, galactic 
diffuse background (gll_iem_v05-revl.fit|14l|) and the 
Second Fermi LAT Catalog (2FGL) sources in ROI of 
10 degrees centered at the position of Table |T] sources. 
The spectrum model was according to the 2FGL. Tar¬ 
get blazars (Table IJ) were fitted with a Log-Parabola 
(LP): 

-(a+/31og(E/Eo)) 

(4) 

because LP is typically used for modeling blazar 
spectrafl^. Where Ab[cm“^ s“^ MeV“^] is nor¬ 
malization parameter, Eg [MeV] is scale parameter, 
— {a -\- /31og(E/Eo)) is spectrum index. If the parame¬ 
ter P is zero, LP is equal to Powar-Law spectrum. In 
this paper we fixed Eg parameter of targets to Eg = 
100 MeV. 

We judged the flux variation of target sources by 
some steps. 

Step I: gamma-ray light curves, which width of the 
time bins was fixed on 30 days (shortest bin in this 
study), was made. 

Step II; if calculated Fbar was not required “Selection 
Criteria”, we adopted more large bin size (60 days, 90 
days, 150 days, and 300 days). 


dN 

'dE 


= Nn 


“Selection Criteria” were as follows, /. More than 40 
% of the calculated integral flux of each bin were de¬ 
tected. II. Significant (over 2cr) variation was detected 
by the test in the analyzed period. If the selection 
criteria I. and II. cleared, Fvar of the target could be 
calculated. 


IV. RESULTS AND DISCUSSION 


Figure [T] shows the light curves of Mrk 421 and PKS 
0454-234 as the light curve samples, which are the typ¬ 
ical HBL and FSRQ sources, respectively. The green 
dashed lines represent the average flux of whole pe¬ 
riod. Each Fvar and averaged flux were culculated 
as Mrk421: Fvar = 39.3 ± 1.3 % averaged flux = 
(2.10 ±0.22) X 10-'^ cm-i s-\ PKS0454: Avar = 58.7 
± 1.3 % averaged flux = (3.07 ± 0.32) x 10“’^ cm“^ 
s“^. Note, the Fvar calculation was performed only 
over the 9 TS bins. 

Another FVarS were obtained in the same method 
and the Fbar as a function of the redshift is shown 
in FigH Square and circle marks indicate TeVCat 
and not TeVCat sources, respectively. Each subclass 
of blazars are plotted in different colors (Blue:FSRQ, 
Red: HBL, Magenta: IBL, Green: LBL, Yellow: FRI, 
Black: uncertain type.). In Figured Fvar indicates 
the variability amplitude of the GeV gamma-ray light 
curve. GeV gamma ray from FSRQs could be detected 
at high redshift(z > 0.5) and have the large Fvar. In 
addition, HBL and IBL assemble in z < 0.5. From 
these features, blazar subclass seems to change along 
the increasing redshift. 

Peculiar features were as follows; 

AO 0235±164 (z = 0.94, LBL) has particularly high 
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FIG. 2: The FVar as a function of the redshift. Square and circle marks indicate TeVCat and not TeVCat sources, 
respectively. Each subclass of blazars are plotted in different colors (Blue:FSRQ, Red: HBL, Magenta: IBL, Green: 
LBL, Yellow: FRI, Black: uncertain type.). 


J^var (117 ± 2.8) in LBLs. This source was discussed 
that it might be FSRQ type blazar 0,0 ; therefore, 
high Fhar value of this source might be caused by the 
FSRQ like characters. S3 0218+357 (z = 0.944, Un¬ 
certain type) which has the highest Fvar (190 ± 2.8) 
in this analyzed sources is a gravitationally lensed 
blazar 03, hence the very high Fvar value might be 
enhanced by the gravitationally lensed effect. 

From Fig. [51 the Fvar as a function of the red- 
shift seems connection with FSRQs —>■ LBLs —>■ HBLs 
(IBLs). This trend shows possibility of the activity 
evolution and blazar class evolution. 

Figure [3] shows the Fhar histogram which is pro¬ 
jected in the vertical axis of Fig. [2] The different 
colors show each subclass of blazars (Blue: FSRQ, 
Red: HBL, Magenta: IBL, Green: LBL, Yellow: FRI, 
Black: uncertain type). 

However, there are some problems in this study. 
First, the middle-high redshift (z > 0.2) low Fvar 
sources were not sufficient for discussions without se¬ 
lection effects. Second, this study could not consid¬ 
ered the short time scale variability (< 30days). Thus, 
it is necessary to analyze the low Fkar sources and 
short time scale. 


V. CONCLUSIONS 

We selected 102 AGNs (blazars) to reveal the rela¬ 
tion between the evolution process of AGNs and the 
blazar sequence. We applied the fractional variability 
amplitude (Avar) to calculate the variability amplitude 
considering the error. 

The analyzed AGNs were selected from the second 
LAT AGN catalog (2LAC)[3 and TeVGat0. The 
analyzed data was Fermi reprocessed pass 7 data be¬ 
tween 2008 August 04 and 2014 June 09, using the 
unbinned likelihood analysis with the Fermi Science 
Tools. 

From these features, blazar subclass seems to 
change along the increasing redshift (connection with 
FSRQs —> LBLs —HBLs (IBLs)). This trend shows 
possibility of the activity evolution and blazar class 
evolution. Peculiar features were as follows; AO 
0235+164 (z = 0.94, LBL) has particularly high Fvar 
(117 ± 2.8) in LBLs. This source was discussed that 
it might be FSRQ type blazar 0 , 0 ; therefore, high 
Avar value of this source might be caused by the FSRQ 
like characters. S3 0218+357 (z = 0.944, Uncertain 
type) which has the highest Avar (190 ± 2.8) in this 
analyzed sources is a gravitationally lensed blazar 0, 
hence the very high Fvar value might be enhanced by 
the gravitationally lensed effect. 
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FIG. 3: The FVar distribution. The width of TVar bin is 20%. The different colors show each subclass of blazars (Blue: 
FSRQ, Red: HBL, Magenta: IBL, Green: LBL, Yellow: FRI, Black: uncertain type.). 


However, there are some problems in this study. 
First, the middle-high redshift (z > 0.2) low Fvar 
sources were not sufficient for discussions without se¬ 
lection effects. Second, this study could not consid¬ 
ered the short time scale variability (< 30days). Thus, 


It is necessary to analyze the low F^ar sources and 
short time scale. 
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